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I Abstract 

We study the free boundary problem for the equations of compressible Euler equations with a vacuum 
boundary condition. Our main goal is to recover in Eulerian coordinates the earlier well-posedness result 
' obtained by Lindblad [TT] for the isentropic Euler equations and extend it to the case of full gas dynamics. 

For technical simplicity we consider the case of an unbounded domain whose boundary has the form of a 
' graph and make short comments about the case of a bounded domain. We prove the local- in-time existence in 

^ ' Sobolev spaces by the technique applied earlier to weakly stable shock waves and characteristic discontinuities 

[51 [5T]. It contains, in particular, the reduction to a fixed domain, using the "good unknown" of Alinhac 
[1], and a suitable Nash-Moser-type iteration scheme. A certain modification of such an approach is caused 
^ i by the fact that the symbol associated to the free surface is not elliptic. This approach is still directly 

applicable to the relativistic version of our problem in the setting of special relativity and we briefly discuss 
\^ . its extension to general relativity. 

(N 



O 

00 

O 

> 



1 Introduction 

Consider the compressible Euler equations with the gravitational field Q G 



><| . dtp + div{pv)^0, (1) 

H ■ 

pi \ dtifv) +div {pvSiv) +\7p = pg, (2) 

dt{p{e+^\v\^))+dW {{p{e+^\v\^)+p)v)^0, (3) 

where p denotes density, w e M"^ fluid velocity, p = p{p, S) pressure, S entropy, and e — e{p, S) internal 
energy. With a state equation of gas, p = p{p,S), and the first principle of thermodynamics, (H])-© is a 
closed system. As the unknown we can fix, for example, the vector U = U{t, x) = (p, v, S). 
We can easily symmetrize system ([I])-© by rewriting it in the nonconservative form 

where — Pp{p, S) is the square of the sound velocity and d/dt = dt + {v, V) (by ( , ) we denote the scalar 
product). Equations ([4]) read as the symmetric quasilinear system 

3 

Ao{U)dtU + ^Aj{U)djU + Q{U) = 0, (5) 



1 



where Q{U) = (0,-^^,0), 
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System ([5]) is symmetric hyperbolic if the the hyperbohcity condition Aq > holds: 

P>0, pp>0. (6) 

One can alternatively consider the isentropic Euler equations, i.e., system ([T]), ([2]) for the same variables 
except for the entropy S. Then, the state equation of gas is p = p{p) and the second inequality in ([6]) is 
understood in the sense that p'{p) > 0. 

We are interested in the motion of an ideal compressible fluid (gas) body in vacuum described by the 
Euler equations (HJ-© (or ([1]), ([2|) for isentropic gas) in a space-time domain n{t) which boundary I](t) = 
{F{t, x) — 0} is to be determined and moves with the velocity of the gas particles at the boundary: 

dF 



dt 



= 0, 



p = 



(7) 



(for all i G [0,T]). This free boundary problem can be used for modeling the motion of the ocean or a star. 
Most results for such kind of problems were earlier obtained for incompressible fluids and the history of 
mathematical studies of incompressible versions of problem (Il])-([3|), ([7]) can be found, for example, in [llj . 

The first result for compressible fluids was obtained by Makino |T3] (see also [H]) who proved the 
local-in-time existence of solutions to problem ((T])-([3]), ([7]) for the case of a polytropic gas and when the 
boundary condition p — Q in ^ is replaced by p = 0. This was done by using a special symmetrization 
of the gas dynamics system that supports vacuum regions. That is, the corresponding symmetric system 
for a new unknown U (see |13l I14| ) is always hyperbolic without assumptions However, employing this 
symmetrization leads to certain non-physical restrictions on the initial data. Therefore, Makino's result 
does not cover the general case. On the other hand, from the physical point of view, the vacuum boundary 
condition pjs = is, of course, more natural than pjs = 0. In particular, ([6]) and ([7]) does not formally allow 
the equation of state of a polytropic gas p =^ ap"' exp(S'/cv). In this connection, as was recommended in 
for the case of boundary condition p\y, — Q one can alternatively think of the pressure as a small constant 
on the boundary (see also Remark 12.11 below) . 

The local-in-time existence for the general case of initial data was recently proved by Lindblad [11] for 
the free boundary problem with non-vanishing density on the boundary for the isentropic Euler equations. 
Namely, the local-in-time existence of smooth solutions of problem 11]), ([2]), fT]) (with ^ = 0) was shown in 
[llj under the natural physical assumption 

dp 
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< -e < 



on 



S(0), 



(8) 



2 



where d/dN = {WF, V), together with the hyperbohcity condition ([6]), provided that the initial domain il{0) 
is difFeomorfic to a baU. The main tool in is the passage to the Lagrangian coordinates for reducing 
the original problem to that in a fixed domain. Such a technique seems most natural for free boundary 
problems with boundary conditions like ([7]). At the same time, for compressible fluids it is connected with 
a lot of technical difficulties and it is not quite clear how to extend the results to similar problems for more 
complicated fluid dynamics models like, for example, relativistic gas dynamics or magnetohydrodynamics. 
Even the extension of the existence theorem in [11] to full gas dynamics does not seem to be just a technical 
matter. 

Remark 1.1 If the domain fl{t) is unbounded, we should additionally assume that the velocity vanishes at 
infinity (as |a;| oo). As follows from the second vector equation in in the absence of gravity {G = 0) 
this contradicts condition ([8]). That is, in the case of an unbounded domain, the presence of gravity is 
absolutely necessary. However, if the domain is bounded, without loss of generality and as was done in [llj . 
the gravity can be neglected as a lower order term (it plays no role in the proof of well-posedness). 

In this paper we propose another approach to studying the well-posedness of problem ([I])-®, ^ (or 
(III), ([2]), 0) and similar free boundary problems for other systems of hyperbolic conservation laws. This 
approach could be probably called "hyperbolic" or "shock waves" approach because it was first applied by 
Blokhin (see [2j and references therein) and Majda [H] to prove the short-time persistence of discontinuous 
shock front solutions to hyperbolic conservation laws. The "hyperbolic" approach to free boundary problems 
does not propose to pass to the Lagrangian coordinates (the more so as this is impossible for shock waves). 
Instead of this we work in the Eulerian coordinates and reduce our free boundary problem to that in a fixed 
domain. More precisely, such a procedure is indeed quite simple if our domain Q{t) is unbounded and its 
boundary has the form of a graph. In this case we reduce our problem to that in a half-space by simple 
straightening of the unknown free surface (for example, a shock front). Otherwise, the technique of reduction 
to a fixed domain is more technically involved (see |12j). but the resulting problem in a fixed domain has no 
principal differences from that for the case of unbounded domains. We can then follow standard arguments 
and reduce the corresponding linearized problem to a linear problem in a half-space by using a fixed partition 
of unity flattering the boundary. Therefore, without loss of generality we can restrict ourself to an unbounded 
initial domain and we do so in this paper. On the other hand, the possibility to treat unbounded domains 
is already a certain advantage of the "hyperbolic" approach. 

Regarding the free boundary problem (HJ-Q, it should be noted that its linearized version is well- 
posed only in a weak sense. It means that the corresponding linear problem satisfies the Kreiss-Lopatinski 
condition but violates the uniform Kreiss-Lopatinski condition [9l [121 IE] ■ This yields losses of derivatives in 
a priori estimates for the linearized problem. Therefore, we are not able to use such estimates to prove the 
existence of solutions to the original nonlinear problem by the fixed-point argument as was done by Blokhin 
or Majda (see also [T^) for uniformly stable shock waves (the uniform Kreiss-Lopatinski condition holds for 
such shocks). Thus, we have to modify the "hyperbolic" approach to apply it to free boundary problems 
whose linearized versions are weakly well-posed. In some sense, this was already done in previous works. We 
should first mention Alinhac's study [1 of rarefaction waves for hyperbolic conservation laws. 

It is well-known that the Nash-Moser method can sometimes compensate the loss of derivatives phe- 
nomenon and to use it we should perform a genuine linearization of our nonlinear problem, i.e., to keep all 
the lower-order terms while linearizing. One of these terms is a first-order term for the perturbation of the 
free surface in the linearized interior equations. To neutralize such a bad term Alinhac proposed to pass to a 
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new unkwnown (so-called "good unknown") and we use this idea for problem (IT])-(l3]), ([7]). Such a technique 
was recently applied to other hyperbolic free boundary value problems. We mean the results of Coulombel 
and Secchi for 2D supersonic vortex sheets and weakly stable shock waves in isentropic gas dynamics 
and author's result for compressible current- vortex sheets [501 [H]. The local-in-time existence of the listed 
weakly stable discontinuities was shown in [5l [21] by a suitable Nash-Moser-type iteration scheme. 

At last, we should note that problem ([T])-([3|), ([7|) is not a quite standard "weakly stable" hyperbolic free 
boundary problem like those studied in [Il[5l[21]. Actually, regardless of the fact that the constant ("frozen") 
coefficients linearized problem for ([IJ-Q, © always satisfies the weak Kreiss-Lopatinski condition, the 
corresponding variable coefficients problem is not unconditionally well-posed and ([8]) is an extra condition 
which is necessary for well-posedness (though, the question on its necessity is a separate and non-trivial 
problem) . This unusual feature is a consequence of the fact that the symbol associated with the free surface 
is not elliptic (see Remark [5^ that leads to a loss of "control on the boundary." Therefore, we have to 
modify somewhat the energy method which we use for deriving a priori estimates for the linearized problem. 
Having in hand a good a priori estimate (so-called tame estimate [1]) for the linearized problem, we prove 
the local existence (and uniqueness) theorem for our nonlinear problem (see Theorem 12.11 below) by the 
Nash-Moser method. 

Such a modified "hyperbolic" approach outlined above allows one to prove a counterpart of Theorem 
12.11 for the relativistic version of problem ([Ij-Q, ([7|) in the setting of special relativity without further 
modifications. Actually, the proof is absolutely the same as for the non-relativistic case and we may drop 
it. Since in the framework of our "hyperbolic" approach we use the energy method (but not the Kreiss 
symmetrizer technique [SJ [T^l [H] ) , the only important point is that the system of relativistic Euler equations 

V„(pu") = 0, V„T"'3 = o (9) 

can be symmetrized (we write down its symmetric form in the last section of the paper). Here Vq is the 
covariant derivative with respect to the metric g with the components gap] p is the particle number density 
in the rest frame (for convenience we use the notations that are consistent with the non-relativistic case); 

h = 1 -\- e ~\- {p/ p) is the specific enthalpy, p is the pressure, e — e{p, S) is the specific internal energy per 
particle, S is the entropy per particle, are components of the four-velocity. The metric g should satisfy 
the Einstein equations. Following [TH] (see also [1|), in the last section of the paper we write down them in 
so-called harmonic coordinates. In the case of special relativity g = diag (—1, 1, 1, 1) and equations ([9]) (in 
the presence of gravity) take the form 

dtipT) + div (pu) ^0, (10) 
dtiphTu) + div {phu (g) u) + Vp = pG, (11) 
dtiphV^ -p) + div {phTu) = 0, (12) 

where 

t:=x°, div:=diva;, x = (x'^ ,x'^ ,x^), u = {u'^,u^,u^), v = {v'^ ,v'^ ,v^) = u/T, 1 + 

r — — {1 — |iip)^^/^ is the Lorentz factor, and the speed of the light is equal to unity. 
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Regarding the free boundary problem for relativistic fluids with a vacuum boundary condition, its local- 
in-time existence was proved by Rendall :18] for the boundary condition — and a special class of initial 
data by generalizing Makino's symmetrization [131 [13 the relativistic case. This result was obtained for 
the setting of general relativity and under the simplifying assumption that the relativistic fluid is isentropic. 
Actually, in the framework of Makino's approach this assumption was just a technical simplification. That 
is, our main goal in this paper is to cover the general case of initial data but for the boundary condition 

As was already noted above, we do not almost need to make efforts for extending Theorem 12.11 to the 
relativistic Euler equations in the setting of special relativity. Concerning the case of general relativity, 
the proof of the existence theorem is based on using harmonic coordinates and the facts that the Einstein 
equations for the metric g can be written in the form of a symmetric hyperbolic system |18j and the metric 
should be smooth on the fluid- vacuum boundary S. More precisely, for the relativistic Euler equations 
we easily obtain a counterpart of Theorem 12. II for any fixed metric, but not only for g = diag (—1, 1, 1, 1). 
Then, roughly speaking, we resolve the relativistic Euler equations by Nash-Moser iterations whereas at each 
Nash-Moser iteration step we find the metric from the Einstein equations by Picard iterations. Actually, we 
do not even need to write down Picard iterations because we know that a unique solution to the Einstein 
equations (for fixed fluid unknowns) written in the form of a symmetric hyperbolic system does exist and this 
is proved by the classical fixed-point argument. Since it makes probably sense to devote a separate paper to 
the case of general relativity we restrict ourself to a schematic proof of the existence theorem. Moreover, we 
do not even formally write down such a theorem in this paper. 

The plan of the rest of the paper is the following. In Section[51 we reduce problem (P)-®, ([7]) to that in a 
fixed domain and state the existence Theorem l2.1l for the reduced problem. In Section[2]we also formulate the 
linearized problem and prove its well-posedness under suitable assumptions on the basic state about which 
we linearize our nonlinear problem ([l])-([3]), ([7]). The main of these assumptions is the physical condition 
([5]) . In Section O for the linearized problem we derive an a priori tame estimate in the Sobolev spaces H'^ 
with s > 3. In Section |4l we first specify compatibility conditions for the initial data and, by constructing 
an approximate solution, reduce our problem to that with zero initial data. Then, we solve the reduced 
problem by a suitable Nash-Moser-type iteration scheme. At last, in Section [5] we describe extensions of the 
result of Theorem 12.11 to special and general relativity. 



2 Basic a priori estimate for the linearized problem 

For technical simplicity (see Remark 12.21 below), we assume that the space-time domain f2(t) is unbounded 
and lies from one side of its free boundary I](i) which has the form of a graph, xi — f(t,x'), x' — (2:2, X3). 
That is, 

VL(t)^{xx> ip{t,x')'\ (13) 

and the function ip{t,x') is to be determined. As for shock waves, using Majda's arguments p2] , we can 
generalize the technique below to the case of an arbitrary compact free surface E. The mapping of il,{t) to 
a fixed domain is just more technically involved when ^l{t) is bounded (see Remark l2.2p . 
For domain (|13p the boundary conditions ([7]) take the form 

dt^ = VN, P = on E(t), (14) 
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and the gravitational field 

g = {G,o,o), 

where vn — {v, N), N — (1, —82^, —d2(p), and G denotes Newton's gravitational constant. Our final goal is 
to find conditions on the initial data 

U{0, x) = Uoix), X e n{0), (^(0, x') = Mx'), e M^ (f 5) 

providing the existence of a smooth solution {U, of the free boundary value problem (fT4|) . (fT5|) in n{t) 
for all t G [0, T], where the time T is small enough. 

To reduce the free boundary value problem ([5]), HH), (fT5|) to that in a fixed domain we straighten, 
as usual, the unknown free surface S. That is, the unknown U being smooth in Q{t) is replaced by the 
vector-function 

U{t,x) := U{t,<^{t,x),x'), 

that is smooth in the fixed domain — {xi > 0, x' Q R^} , where $(t, 0,a:') — (p{t,x') and 9i$ > 0. As 
in |2f| . to avoid assumptions about compact support of the initial data in the nonlinear existence theorem 
and work globally in M.^ we use the choice of $(t, x) similar to that suggested by Metivier [15] : 

x) := xi + x), x) xi^iMt, a;'), 

where x G equals to 1 on [0,1], and ||x'||loc(R) l/-^- Then, the fulfillment of the requirement 

9i$ > is guaranteed if we consider solutions for which ||</'||l^([o,t]xr2) 5; 1- The last is fulfilled if, without 
loss of generality, we consider the initial data satisfying ||</5o||loo(R^) — l/^j ^-nd the time T in our existence 
theorem is sufficiently small. 

Dropping for convenience tildes in U, we reduce ([5]), (O, (US]) to the initial boundary value problem 

L(C/,*) = in[0,T]xRi^, (16) 

B(C/, (/?) = on [0, T] X {xi = 0} x M^, (17) 
U\t=o = Uo mRl, ^\t=o = Vo inK', (18) 
where L(C/, *) = L{U, «-)[/ + Q{U), 

L{u, *) = Ao{u)dt + Ii(c/, *)ai + A2{u)d2 + A3(c/)a3, 
^i(c/,*) = ^^(^le^) - MU)dt^ ~ ^ Afe(c/)afeVi/) 

^ fc=2 

(9i<i> = 1 + 9i\l/), and (fT7)) is the compact form of the boundary conditions 

dt^-VN = p = on [0, T] X {xi = 0} X R^. 

We are now in a position to state the local- in-time existence theorem for problem (fT6|) - p8)) . Clearly, this 
theorem implies a corresponding theorem for the original problem ^ , (jl4p , (|15p . 

Theorem 2.1 Lei to € N anrf to > 6. Suppose the initial data \10(l . with 

([/o-?7,^o)eH™+'(Ki)xi/"+^(M2) and p(po, ^0) - ei G H™+7(M^), 
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satisfy the hyperbolicity condition ([6]) for all x G M'^ and are compatible up to order m + 7 in the sense of 
Definition \4.1\ Here 

(26X1,0,0,0,0), ei=2e/G (e = const > 0). 
Let also the initial data satisfy the physical condition 

dip>e>0 at xi = (19) 

for all x' G M^. Then, there exists a sufficiently short time T > such that problem hlb]) - lT^) has a unique 
solution 

{U,(p) G {t/ + i7™([0,T] X Rl)} X H"\[0,T] x R^). 
Moreover, p - ei E iJ™([0,T] x Rl). 

Remark 2.1 The hyperbolicity condition ([6]) which should be satisfied for all x G R\_ implies that the 
function po{x) = p{po, So){x) cannot vanish at infinity. Indeed, in Theorem 12.11 we assume that po — ei G 
ff™+^(Mi^). On the other hand, ([6]) together with the boundary condition p\-s = do not formally allow 
the equation of state of a polytropic gas (or a 7-law gas for isentropic gas dynamics). However, as was 
noted in [llj . from a physical point of view we can alternatively think of the pressure as a small positive 
constant e on the boundary. One can easily generalize the result of Theorem 12. II to the case of the boundary 
condition p\xi=o — £• More precisely, we now assume that Uq — U — Cq G i?™+^(Mi^) and prove that 
U -U -Cq e H"'{[0, T] X R\), where Co = (e, 0, 0, 0, 0). Indeed, making the change of unknown p' =p-e 
and omitting the primes, we obtain problem (|16p - p^ with the matrices Aa{U + Co). The further arguments 
are almost the same as in the proof of Theorem 12.11 (see below). 

Remark 2.2 Inequality is a counterpart of the physical condition ([5]) for the unbounded domain (fT5|) . 
If the domain is bounded and its initial boundary S(0) is a compact co-dimension one surface in M^, we can 
follow Majda's arguments [12] (see also [21 sect. 12.4.2]). More precisely, we can make {locally in time) a 
change of variables that sends all boundary locations E(t) to the initial surface 5](0). We refer the reader to 
[T^ [5] for details of such a change of variables. In particular, it requires the application of the Weingarten 
map while writing down boundary conditions on S(0). The resulting initial boundary value problem is a 
problem in the fixed domain $1(0). Its principal difference from problem (ll6 |) - ([T8l) is that we have to deal 
with a problem in a fixed compact domain instead of a half-space. For this problem the proof of a counterpart 
of Theorem 12. H is more technical, but the ideas are basically the same as for Theorem 12. II For instance, we 
should reduce the corresponding linearized problem to that in a half-space by using a fixed partition of unity 
flattering the boundary. The resulting linearized problem in a half-space will not have principal differences 
from the linearized problem for (fT6 l) -([T8 l) . Only its coefficients will be more technically complicated than 
those for the linearization of (|16 p -p ^ . Therefore, as is usually done for shock waves or other types or strong 
discontinuities (see, e.g., [HOdl]), in this paper we restrict ourself to the case of an imbounded domain 
whose boundary has a form of a graph. 

The existence of solutions in Theorem 12.11 will be proved by Nash-Moser iterations. The main tool for 
proving the convergence of the Nash-Moser iteration scheme is a so-called tame estimate [1] [5l [21] for the 
linearized problem. In this section, we derive a basic a priori L2^cstimatc for the linearized problem by 
the energy method. This estimate is a basis for deriving the tame estimate in Sobolev spaces (see the next 
section) and implies uniqueness of a solution to the nonlinear problem (|16p - (|18p that can be proved by 
standard argument. 
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Let us first pass to the new unknown U' = {p',v,S) = U — U. For U' system (fT6)) is rewritten as 

h'{U', -f) := L{U' + U, *)C/' + A^(U' + U, 'f)diU + Q{U' + f/) = 0, 

where diU = (2e, 0, 0, 0, 0). Let p'{p',S) := p(p,S'), A'^{U') := Aa{U), Q'{U') := Q{U), and := Uq-^U. 
Then, omitting the primes, for the new unknown we get the system 



UU, *) := L{U, *)[/ + A^{U, + Q{U) = in [0, T] x 



(20) 



with the boundary conditions ([T7]) and the initial data (1181) . From now on we will work with problem 
(1^ . (fTT]) . HI]). We should now prove the existence of its solution, U e iJ™([0,T] x R^), assuming that 
Uq E iJ™+^(Ri^). For the initial data for the new unknown we assume that 



This guarantees the fulfillment of assumption ()19|) for the original unknown. 
Remark 2.3 We easily compute the boundary matrix: 



(21) 





( 


_L 


1 


-92* 


-93* 
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/ 



where f — vi — W2(?2* — fs^s* — 9t*. The vector-function Ai{U, *°)9iC/ cannot belong to a Sobolev space 
on M.^ because its second component is 2e/(9i$). However, if problem ([20]) . (fT7|) . (fT8|) has a solution from a 
Sobolev space and Theorem 12 . II takes place, then the sum Ai{U, *)9i?7 + Q{U) already belongs to a Sobolev 
space because 

^-Gp^ -G{p 6i) - 2e|| e H"\[0, T] x ). 
Thus, for our case of an unbounded domain the presence of gravity is of great importance (see also Remark 

ini). 

We now formulate the linearized problem. Consider 



-oo,T] 



dilr := (-oo,T] x {xi = 0} x 



Let 



be a given sufficiently smooth vector-function, with U ~ {p,v,S), and 



(22) 



(23) 



where ii' > is a constant. Moreover, without loss of generality we assume that ||<p|| L<^(anT) ^ 1- This 
implies 9i$ > 1/2, with <I>(t, a;) :— xi +*(t, x), *(i, x) := x{^i)'f{tj x'). We also assume that the basic state 
about which we shall linearize problem (PO)) . p?)) satisfies the hyperbolicity condition ^ in Q.t, 



p{p,S)>0, pp{p,S)>0, 



(24) 
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the first boundary condition in (jl7p . 
and the assumption (pij) . 



dtifj - vn\xi=o = 0, 
9ipUi=o > -e, 



(25) 
(26) 



where vn — vi — V2d2'P — v^d^ip. 

The linearized equations for (|20p and (|17p for determining small perturbations {5U, Sip) read (below we 
drop S): 

L'(C/, $)([/, L(C/, $)[/ + CiU, - {L(C/, $)*}M^^±^ ^ 

+ V2d2(p + vsdaifi ~ vn \ 
I 

where vpf = vi ~ V2d2'p — v^d^ip, and the matrix C{U, 5*) is determined as follows: 



CiU, *){/ = ((7, VuAo{U))dtU + (U, VMU, ^))diU 



k=2 



( \ 

-5Pp(p, S')p - gps(p, 5)5 






(dA{Y) 



Y = (yi, . . . ,y5)- 



Here, as usual, we introduce the source terms f — (/i, . . . , /a) and g — ((?i, 52) to make the interior equations 
and the boundary conditions inhomogeneous. 

The differential operator L'(C/, ^) is a first order operator in 5* = x{xi)'p{t,x'). Following Alinhac [T] 
and introducing the "good unknown" 



^di{U + U), 

9i$ 



(27) 



we simplify the linearized interior equations: 



L{U, *)t/ + C{U, *)[/ di{h{U, ^)} = /. 

As in [U El [20I [21] , we drop the zero-order term in 5" in ([28ll and consider the effective linear operators 
V^{U, $)J7 := L(t/, $)J7 + C(C/, $)f7 

= Ao{U)dtlj + MU, d)diU + A2{U)d2U + A:i{U)d^U + C{U , $)J7 



(28) 



(29) 



In the subsequent nonlinear analysis the dropped term in l|28p will be considered as an error term at each 
Nash-Moser iteration step. 

Regarding the boundary differential operator B', in terms of unknown ((27|) it reads: 



?',([/, ^)({7,(^) :=B'(&,^)(C/,^) = 



■p + ip{2e + dip) 



(30) 
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where vyi ^ vi — V2d2<p — vsda'fi- Thus, the hnear problem for (U , ip) has the form 

{U, (p)^0 for t < 0, 



(31) 
(32) 
(33) 



where / and g vanish in the past. We consider the case of zero initial data, that is usual assumption, and 
postpone the case of nonzero initial data to the nonlinear analysis (construction of a so-called approximate 
solution). 

On the basic state the boundary matrix Ai has the form 



Ai(J7,4') 



where 



In view of (gS]), 



p(p, S), 





( ± 

pc? 


1 




C»3$ 






1 


1 


/5f 








































/5f 




















f ) 




= Pp{p,S), f 


= Vl 


- ^292* 










= Vn\xi-- 


=0 — 


dtv = 0. 







We see that the boundary matrix Ai{U , 5*) is singular on the boundary xi ^ Q (it is of constant rank 2 at 
Xi = 0). That is, (|3ip - (|33p is a hyperbolic problem with characteristic boundary of constant multiplicity. 

It is convenient to separate "characteristic" and "noncharacteristic" unknowns. For this purpose we 
introduce the new unknown 

V = {p,Vn,V2,V3,S), 

where v„ = vi - V2d2^ - vsda'i' {vn\xi=o = vn\xi=o)- We have [/ = JV, with 

/ 1 \ 

1 ^2$ ^3$ 

1 

1 

1 



J 



\ 



J 



Then, system pip is equivalently rewritten as 

3 

Ao{U, ^)dtV + MU, ^)dkV + AiiU, $)F = T{U, $) 



(34) 



fc=i 



where Aa — AaJ {a ~ 0,2,3), Ai = J~^AiJ, T = f ■ The boundary matrix Ai in system ([M]) 
has the form 



Ai = —^A(i) + A(i^), A(i) 



/ 





1 











\ 




1 






















































V 

















/ 



^(o)Ui=o = 0, 



(35) 
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i.e., Vn = {PtVu) is the "noncharacteristic" part of the vector V . The exphcit form of A[o) is of no interest, 
and it is only important that, in view (|25p . -4(o)|2;i=o ~ 0. The boundary matrix Ai on the boundary a;i = 
has one positive ("outgoing") eigenvalue. Since one of the boundary conditions is needed for determining 
the function ip, the correct number of boundary conditions is two (that is the case in (1321) ). Hence, the 
hyperbolic problem (|3ip - p3|) has the property of maximality |16j . 

By standard argument we get for system ([3T|) the energy inequality 

I{t)~2j^^ pm\.,=odx'ds<C{K)(^\f\\l^^,,^)+ /(s)ds^, (36) 

where I{t) = J^s {AoV, V) dx and C — C{K) > is a constant depending on K (see p3l)). In view of the 
boundary conditions ([32|l . one has 

-2pvn\xi=o = 2(^30 - g2)idt(p + V2d2(p + hdsV ~ 'pdiVN ~ 5i)Ui=o 
= dt {a\x^=oV^ - 2.g2¥'} - {dto. + d2{v2a) +93(^30) - 2ddiVN}\xi=o 
+2 {dtg2 + d2{v292) + dz(v'i92) + g2diVN - gia} \xi=o V + 2.9ig2 
+82 {v2cnp^ - 2v2g2ip} + ds {vsaip'^ - 2v3g2(p} , 
where a = 2e + dip. Then, using the Young inequality, from (|36[) we obtain 

m + J^^ (2e + dip\x,=o) dx' < C(i^){||/||i^(j,^) + \\g\\h(an^) 

+ J* (/(s) + ||^(s)||i^(R.)) ds}. 

Taking into account assumption (j26p and applying Gronwall's lemma, we finally deduce the basic a priori 
L2-estimate 

\\U\\L2(nT) + \\'P\\L2idnT) < C{K) {\\f\\L2(nT) + 11.911^1(00^)} • (37) 

Remark 2.4 In the a priori estimate ([57|) we have a loss of one derivative from the source term g to the 
solution (more precisely, we loose one derivative only from 172 but not from gi). This is quite natural 
because one can check that the constant coefficients linearized problem, i.e., problem (I3ip -(l33 |) with frozen 
coefficients satisfies the Kreiss-Lopatinski condition but violates the uniform Kreiss-Lopatinski condition 
[HI [H] • Although the weak Kreiss-Lopatinski condition holds we had to assume the fulfillment of the extra 
condition (|26p while deriving the a priori estimate (|37p . This is very unusual for hyperbolic initial boundary 
value problems because, as a rule (see, e.g., [5J [3T]), the fulfillment of the Kreiss-Lopatinski condition is 
enough for obtaining a priori estimates. Actually, in our case the appearance of an extra condition on the 
level of variable coefficients linear analysis is caused by the fact that the symbol associated to the free surface 
is not elliptic, i.e., we are not able to resolve our boundary conditions (I32|) for the gradient {dfif, d2ip, d^tp). 
Therefore, it is also natural that in estimate ([37|) we "lose one derivative from the front", i.e., we do not 
have the H^-norm of (p in the left-hand side of ((37|l . 

Since in estimate ([37]) we do not lose derivatives from the source term / to the solution, the existence of 
solutions to problem (|3i p - (p3|) can be proved by the classical argument of Lax and Phillips [10]. Indeed, we 
first reduce our problem to one with homogeneous boundary conditions by subtracting from the solution a 
more regular function (see, e.g., lUj). Namely, there exists U — {p,v,S) G i?''+^(riT) vanishing in the past 
such that 

-VN = gi, P = g2 ondflr, 
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where v-^ — vi — V2d20 — v^d^ip. liU = U'^ + U , then satisfies ((3T|) ~([33 | with .g = and / ~ f \ where 
= / — Lg ([/,'!')[/. That is, it is enough to prove the existence of a solution {U,ip) to problem ((3T|) -([33 | 
with (7 = 0. For this problem we have the estimate 

||f>l|L.(n.) + ML.ianr) < C{K) ||/|U,(n.). (38) 

Having in hand estimate (|38p with no loss of derivatives we may use the classical argument in |10j . In 
particular, we define a dual problem for ([3T|l -([33 l) as follows: 

K{U,^YU = f inrir, (39) 
dtp + d2{v2p) + ds^vsp) + pdiVN + vno- = on dflr, (40) 
[7 = for i < 0, (41) 

where U — (p, v, S), = vi — V2d2'f — v^dy,(p, and 

K* = -K + C + - dtAo - diM - d2A2 - dsAs. 

Problem (l39|)-(l41]) is indeed a dual problem for (l31])-(l33l) because for aU U £ F^^t) and U G H'^i^r), 
with U\t=T = 0, satisfying the homogeneous boundary conditions (f32|) (with g = 0) and (|40|) respectively, 
one has 

(L^[/,[7)i,(a,) - (;7,L^*C7)i,(o,) = -(Ii?7,t7)i,(go,) = -(AT^,^)l.(<50,) = 0, 
where V = J~^U. For the dual problem ([39)) " (|4T|) we can easily get the inequality 

m + 1 ^^^^/U^o<^^' < ciK) [mu,^,+J^ (lis) + iip-,...„(.)iii.(«.,) d.} 

{I{t) = /jj3 {AaV, V) dx) which, in view of condition ((26l) . implies the L2-estimatc 

We omit further arguments which are really classical and refer to [TU] (see also, e.g., [Hill]). Thus, we have 
the following well-posedness theorem for the linearized problem ([3T|) - p3|) . 

Theorem 2.2 Let assumptions P5)) -(|26p are fulfilled for the basic state i2^) . Then for all (/, g) G L2(f^T) X 
H^{dflT) that vanish in the past problem ll31]) - ^33\) has a unique solution [U , ip) G L2{flT) x L2{dflT)- This 
solution obeys the a priori estimate (|37p . 

Remark 2.5 Strictly speaking, the uniqueness of the solution to problem (|3ip - (p3)) follows from estimate 
([57)l . provided that our solution belongs to _ff^(riy) x H^{dftT)- We omit here a formal proof of the 
existence of solutions having an arbitrary degree of smoothness, and we shall suppose that the existence 
result of Theorem 12. 21 is also valid for the function spaces H^{flT) x H''{dilT), with s > 1. In this case exact 
assumptions about the regularity of the basic state will be made in Sect. 3, where we prove a tame a priori 
estimate in iJ^(f2T) x iJ*(9f2T) with s large enough. 

3 Tame estimate for the linearized problem 

We are going to derive a tame a priori estimate in H" for problem (PT|) - ([M)) . with s large enough. This tame 
estimate (see Theorem 13.11 below) being, roughly speaking, linear in high norms (that are multiplied by low 
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norms) is with no loss of derivatives from /, with the loss of one derivative from and with a fixed loss 
of derivatives with respect to the coefficients, i.e., with respect to the basic state (f22|) . Although problem 
pip - (|33|) is a hyperbolic problem with characteristic boundary that implies a natural loss of control on 
derivatives in the normal direction we manage to compensate this loss and derive higher order estimates in 
usual Sobolev spaces. This is achieved by using the same idea as in [191 E] and estimating missing normal 
derivatives through a vorticity-type linearized system. 

Theorem 3.1 Let T > and s e N, with s > 3. Assume that the basic state {U,ip) S H^^^{VIt) x 
H^^^{d^T) satisfies assumptions (23 \) -^26 \) and 

\\U\\minT) + \MHHdnT) <K, (42) 

where K > is a constant. Let also the data {f,g) G H'^{Q,t) x H'^^^{dQ.T) vanish in the past. Then there 
exists a positive constant Kq that does not depend on s and T and there exists a constant C{Kq) > such 
that, if K < Kq, then there exists a unique solution (U,ip) € H''(JIt) x H^(dil,T) to problem I131 \} - I[33\) that 
obeys the a priori tame estimate 

\\U\\H-{nT) + MH-idnr) ^ C;iKo)\ \\f\\H-{nT) + WgWH^+^dnr) 

(43) 

+ {\\f\\H^nT) + \\fj\\HHdnT)){\\U\\H-+^(nT) + W^PWH^+^idnT)) \ 
for a sufficiently short time T . 

Proof. Since arguments below are quite standard we somewhere will drop detailed calculations. By applying 
to system (|34|) the operator d^^^^ = 9"° 9^^, with \a\ = |(q;o, Q!2, aa)] < s, one gets 

/ ( A^t^an^, d^^^V)dx - 2 / 9ra„p 9ra„«jv U,=o Ax' ds = 7^, (44) 

where 

n = / ({divA9,^„y - 2^[9r,„, - 2d':,MiV) + 2d^,^T],d:,,,v)dxAs, 

divA = J2'j=o'^j^j (^0 ■— dt), and we use the notation of commutator: [a,b]c :— a{bc) — b{ac). Using the 
Moser-type calculus inequalities 

< C {\\u\\H^{nT)MLU<^T) + \\u\\LU<^T)\MH^{nT)) > (45) 

\\F{u)\\H^^n^^ < C(M) (1 + MHsi^n^)) , (46) 

where the function F is a C°° function of u, and M is such a positive constant that ||m||l^(Ot) — -^j '^^ 
estimate the right-hand side in ([44]): 

n<c{K)[\\v\\ 

ll/llL(o.))(l + l|cocff 11^+0 (47) 

with |lcoeff|l„ := \\U\\H'^{nT) + ll<^llH-(9f^T)- 

Taking into account the boundary conditions, we have: 

-^te.nP^te.n'"N\x^=0 = '^dZ.ni^a - g2)d^^^{dtip + V2d2V + V^d^ip -ipdiVN " 5l)Ui=0 
= dt {aUi=o {dte^nVf - 29ran5'2 d^^^ifi} +...+ 82 {V2\x^=i)[dte,n, ^lPUi=o]y'} d^^^Lp + . . . , 
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where the underhned term is just a typical one that gives a biggest loss of derivatives from the coefficients 
in the final a priori estimate (|43|) . Indeed, using the calculus inequality (|45p and the trace theorem, we get 

||52{t)2Ui=0[9t"an.5lPUl=0]'^}|lL(R2) < W { 1 1 1 1 ff. (go,) + 1 1 <^ 1 1 1^ (SJl^) (l + 1 1 Ui =0 1 1 ^. + 2 (g^^) ) } 

Omitting detailed calculations, from ((44l) and (|47)) we obtain 

III^WIIIL,. + III^WII?/=(M2) < C{K)Mit), (48) 
where ^ 



^{T) = + ll5ll?^.+i(9o.) + + + ll/lli^(o.)) (1 + llcoeffll^^s) , 

lll"WIII?an,™ := E l|5ran"(t)llL(M3^), II I" W I II (D) := E • • ^^(*) 1 1 (15) 
Q|<m i=0 

(Z) = M? or I? = M"^). Since only the biggest loss of derivatives from the coefficients will play the role for 
obtaining the final tame estimate, we have roughened inequality pS)) by choosing the biggest loss. 
It follows from ([M]) and that 

3 

0, 0, 0) = (ai$)^(i) (t - AodtV ~ E ^^^^"^ - -^4^^ - A(o)diV) . (49) 
Applying to (I49p the operator d^^^, with |/3| < s — 1, using decompositions like 

dlUBd^v) = Bdt,Av + [oL, Bm, 

taking into account the fact that ^(o)|ki=o = 0, and employing counterparts of the calculus inequalities (I45p 
and ((46)) for the "layerwise" norms |||(-)(^)||| (see [19]), one gets 

ii9i5fa„K(<)iii ) < ciK){\\\vm\i,,^ + \w^^^Lvml,^s ) + iwvmi.-r^u^ . 

+ III/(^)III^»-i(r3 ) + yWWwi^inT) + ll/llL(f..)) (1 + Il|coeff(t)|||?^0}^ 

where a = (7{xi) G C°°(M+) is a monotone increasing function such that (t{xi) = xi in a neighborhood 
of the origin and cr(xi) = 1 for xi large enough. Since (7\xi=o = we do not need to use the boundary 
conditions to estimate cr9^ 97an^i with j + I7I < s, and we easily get the inequality 

+ + ll/llL(n.)) (1 + l|coeff||?+0}. 

Taking into account Sobolev's embedding in one space dimension, 

lll'"(*)lllff'"-i(D) ^ ll"llLoe([0,t],ff™-i(Ci)) - C'||u||^™([o,t]xD)i 

and combining p8)) . ((50)) . and ((5T|) for j = 1, we obtain 

fc 

i?^MK=)+E E ii5i9ra„K(i)iiL(K3^)<c(x)Ai(t), (52) 



(51) 



2 

tan,s 



= 1 \a\<s-i 
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with fc = 1. 

Estimate (|52|) for fc = s is easily proved by finite induction and equivalently rewritten as 



2 

tan,s 



lK(i)||^e(R3 ) + < C{K)M{t). (53) 



Missing normal derivatives in for the "characteristic" part ('i'2, ^3, S) of the unknown V can be estimated 
from the last equation in pT|) . 

dtS+^Uw,V)S + {u,\/)s\^ h, (54) 
and a system for the linearized vorticity ^ = V x 5, where 

V ^ {vi,Vr2,VT3), Vrk^{v,Tk), T2 = (a2(^, 1, 0), Ta = (93(^, 0, 1), 

This system is obtained by applying the curl operator to the equation for v following from pip . 

9i$ ( p{p,S) J 

{fv = (/2, /r2 , /ra), /r^ = {fv,Tk), fv = (/2, /s, /4)), and has the form 

Ct + ^(«),V)e + l.o.t = V X /„, (55) 

where l.o.t. are lower-order terms which exact form has no meaning. 

Both equations ([51)) and (|55p do not need boundary conditions because, in view of (|^5|) . the first compo- 
nent of the vector w is zero on the boundary xi = 0. Therefore, omitting detailed calculations and combining 
corresponding estimates for the normal derivatives of the "characteristic" unknown (i'2,1'3, S) with (|53p . we 
deduce the inequality 

I{t) < C{K) |aA(T) + ^ I{s) dsj . 
Applying then Gronwall's lemma, one gets 

X{t) < C{K) e^(^)^AA(T) 
(X(0) — 0, see ([55]) ). Integrating the last inequality over the interval [0,T], we come to the estimate 

WrH.i^n.) + \WrH^(on.) < C(X)Te^W^^(r). (56) 

Recall that U — JV. Taking into account the decomposition J{(p) = I + Joi0) and Jo(0) = 0, using (|^ 
together with the improved calculus inequality (|46|) for the case F{0) — 0, 

\\F{u)\\Hs^nr)<CiM)\\u\\Hs^nr), 
and applying Sobolev's embedding in one space dimension, we obtain 

mh^n.^ - WV + JoVWl.^^^^ < C{K){\\V\\l.^^^^ + rilL(o.)l|coeff||?) 

(57) 
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Inequalities (|56|) and ((57)l imply 

ril?^=(o.) + ll^lll^M^o.) < C(if)Te^W^AA(r). (58) 

Taking into account Theorem 12.21 and Remark |2.5[ we have the well-posedness of problem (PT|) - ([55)l in 
H''{nT) X H''{dilT)- Applying Sobolev's embeddings, from (|58|) with s > 3 we get 

+ {\\U\\H^nT) + WfWH^dnT) + WfWH^nT)) (l|f^ll-f/=+3(nT) + ll<Pll-ff=+3(aor)) p 

where we have absorbed some norms ||C/||if3(Q^) and |l</3||ij3(gQ^) in the left-hand side by choosing T small 
enough. Considering ([59| for s = 3 and using (|42p. we obtain for T small enough that 

\\U\\H-^nT) + MH^ianr) < C{Ko) {||/||h3(o^) + WgWrnianT)} ■ (60) 

It is natural to assume that T < 1 and, hence, we can suppose that the constant C'{Kq) does not depend on 
T. Inequalities ^ and ^ imply □ 



4 Nash-Moser iteration 

To use the tame estimate (|43)) for the proof of convergence of the Nash-Moser iteration, we should reduce 
our nonlinear problem (PO)) . p?]) . (fTS]) on [0,T] x Mj^ to that on fix which solutions vanish in the past. 
This is achieved by the classical argument suggesting to absorb the initial data into the interior equations 
by constructing a so-called approximate solution. Before constructing the approximate solution we have to 
define compatibility conditions for the initial data (jlSp. 

(Uoyifo) = (PO, "1,0, 1'2, 0,^3,0, -S'c'Po)- 

Assuming that the hyperbolicity condition ^ is satisfied, we rewrite system in the form 

dtU = ~{Aa{U))-^ (li(C/, •^)diU + A2{U)d2U + A3{U)d3U + A,{U, ^)diU + Q(C/)] . (61) 
The traces 

Uj = {pj,vij,V2,j,V3j,Sj) = d'tU\t=o and = 9tV|t=o, 

with J > 1, are recursively defined by the formal application of the differential operator d{~^ to the boundary 
condition 

dt(p = {vi - V2d2(p - v3d3(p)\^,=o (62) 
and (|6ip and evaluating d^ip and dfU at i = 0. Moreover, '^j = 9j5'|t=o — x{^i)'Pj- 

We naturally define the zero-order compatibility condition as poUi=o — 0- Note that, unlike the case 
when the symbol associated with the free surface is elliptic [3 [151 EI] , this condition does not contain the 
function ipQ. Evaluating ([62]) at t = 0, we get 

^1 = (wi,0 - V2Md2V>0 - V3fid3ipo) \xi=0, (63) 

and then, with dt^\t=o '■= = xi^i)^^!^ from evaluated at i = we define Ui. The first-order 
compatibility condition pi\xi=o = will implicitly depend on ipo and ipi. Knowing ipi and Ui we can then 
find (/?2, U2, etc. The following lemma is the analogue of Lemma 4.2.1 in [15] . Lemma 2 in [3], and Lemma 

5 in [21]. 
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Lemma 4.1 Let E H, ^ > 3, Uq £ H'^{M.^), and (pa E i7''(R^). Then, the procedure described above 
determines Uj G H^^^^M.\) and ipj G iJ^^-' (R^) for j = 1, . . . , fi. Moreover, 

{\\Uj\\H>^-J(Rl) + WfjWH^'-HM^)) < CMq, (64) 

i=i 

where 

Mo = ||C^oll_fi-M(R^) + il</'o||fft'(R2), (65) 
the constant C > depends only on /i and the norms \\Uo\\w'^ {R^) '^'^'^ ll'PollvFi (R^)- 

The proof is almost evident and based on the multipUcative properties of Sobolev spaces (Remark 12.31 
should be also taken into account). 

Definition 4.1 Let ^ e N, /i > 3. The initial data (Uq, tpo) E iJ^(Mij_) x iJ''(R^) are said to be compatible 
up to order when {Uj,(pj) satisfy 

PjUi=o = (66) 

for j = 0,..., fi. 

We are now ready to construct the approximate solution. 

Lemma 4.2 Suppose the initial data il8\) are compatible up to order /i and satisfy the assumptions of 
Theorem [2J\ (i.e., for all x eR\ and l[2l\) ). Then there exists a vector-function ([/", ip"-) E i/^+^(ilT) x 
H^^^{dilT), that is further called the approximate solution to problem (|20p . p7|) . ()18p . such that 

diUU%^'')\t=o^Q forj^Q,...,fi-l, (67) 

and it satisfies the boundary conditions jlT^ , where = xi^i)^"' ■ Moreover, the approximate solution 
obeys the estimate 

||C/1Ih.+i(ji,) + < Ci(Mo) (68) 

and satisfies the hyperbolicity condition (0) on VIt as well as condition i21\) on d^T , where Ci = Ci{Mq) > 
is a constant depending on Mq (see H65\) ). Moreover, p°- — ti — p{p°',S°') — ei G H^^'^{Q.t). 

Proof. Consider functions U" E i/^+i(M x and Lp" E H>'+^{M?) such that 

dl C/" I t=o = Uj E H^-^ {Rl ) , di ^'^ I t=Q = 'Pj E H^^-^ (R2 ) for J = 0, . . . , M, 

where Uj and fj are given by Lemma HTTl Thanks to (pS]) and ([55]) we can choose C/° and ip'^ that satisfy the 
boundary conditions (fT7| . By using a cut-off function we can suppose that (J7", (/9°) vanishes outside of 
the interval [-T,T], i.e., {V^^p'') E Hi^+'^{nT) x Ht^+^idflT)- Applying Sobolev's embeddings, we rewrite 
estimate ([64|) as 

{WUjWh^-HrD + Mn^-Hm) < ^(^^o), (69) 

where C = C{Mq) > is a constant depending on Mq. The estimate (|68p follows from (|69|) and the 
continuity of the lifting operators from the hyperplane t = to R x R'^ . Conditions (|67p hold thanks to the 
properties of {Uj,(pj) given by Lemma [4.11 At last, since {U°',tp"') satisfies the hyperbolicity condition ([6]) 
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and condition (|2ip at < = 0, in the above procedure we can choose ([/", </?°) that it satisfies ([6]) and l|2ip for 
all times t e [~T,T]. The condition — ci e i/^+^(fiT) is just an assumption on the state equation. □ 

Without loss of generality we can suppose that 

I|C^o||hm(r3^) + ||</5o||hm(k2) < 1, IIvjoHhm (R2) < 1/2. (70) 

Then for a sufficiently short time interval [0, T] the smooth solution which existence we are going to prove 
satisfies ||<y5||Lo„([o.T]xR2) < 1 that implies 9i$ > 1/2 (recall that ||x'||l^(r) < 1/2, sec Scction[2]). Let /i is 
an integer number that will appear in the regularity assumption for the initial data in the existence theorem 
for problem ([20|) . p7|) . (fT8|) . Rmming ahead, we take /i = rn + 7, with m > 6 (see Theorem 12. ip . In the end 
of this section we will sec that this choice is suitable. Taking into account ([70]) . we rewrite (pS]) as 



where C* = Ci(l). 
Let us introduce 



l|t/1k™+«(o.) + ll^1li/"-+s(an.) < a, (71) 



-LJ/'^,*'') fort>0, , , 

I fort<0. 



Since (C/",(p'') e H''''+^{VLt) x H'^+^idnr), using dST]), we get /"^ e i7"+^(f7T) and 

|iril//".+^(o,) <'5o(T), (73) 

where the constant (5o(T) — > as T 0. The crucial role in the proof of the fact that /" belongs to a Sobolev 
space is played by the presence of gravity (see Remark l2.3p . To prove estimate ([73| we use the Moser-type 
and embedding inequalities and the fact that /° vanishes in the past. Then, given the approximate solution 
defined in Lemma l4?2l {U,(p) = {U°-,ip°') + {U,'f>) is a solution of the original problem (l20l) . (fT7| . (fT8| on 
[0, T] X Ri^ if ([/, (^) satisfies the following problem on ^It (tildes are dropped): 

C[U,^)^r inf^T, (74) 
B{U, (^) = on OVLt, (75) 
(C/, vs) = for t < 0, (76) 

where £(1/,^') := L(C/'^ + [/, ^-"^ + 5-) - L([/'^,*'^), B{U,ip) := B([/" + [/, (^'^ + <^). From now on we 
concentrate on the proof of the existence of solutions to problem (|74p - ([75|) . 

We solve problem (|74p -(l76 l) by a suitable Nash-Moser-type iteration scheme. In short, this scheme 
is a modified Newton's scheme and at each Nash-Moser iteration step we smooth the coefficient u„ of a 
corresponding linear problem for 5un = u„+i — ii„. Errors of a classical Nash-Moser iteration are the 
"quadratic" error of Newton's scheme and the "substitution" error caused by the application of smoothing 
operators Se (see, e.g., [7J and references therein). As in [3 [21], in our case the Nash-Moser procedure is 
not completely standard and we have the additional error caused by the introduction of an intermediate (or 
modified) state m„+i/2 satisfying some nonlinear constraints. In our case, the main constraint is condition 
(1^ that was required to be fulfilled for the basic state (j^ . Also the additional error is caused by dropping 
the zero-order term in ^E* in the linearized interior equations written in terms of the "good unknown" (see 
(|27|) - ((29l) ). We first hst the important properties of smoothing operators [Tl[5l[7]. 
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Proposition 4.1 There exists such a family {S'ejeoi of smoothing operators in H''{Qt) acting on the class 
of functions vanishing in the past that 

\\Seu\\H^^nr)<Ce^^'"'>^\HH^inT), a,/3>0, (77) 
\\Sgu-u\\Hf^^n^^<Ce^-''\\u\\H^^nr), 0</3<a, (78) 

\\-^Seu\\^,^^^^<Ce^'-"-'\\u\\H^^n^-^, a,(3>0, (79) 

where C > is a constant, and (/3 — a)+ := niax(0, P — a). Moreover, there is another family of smoothing 
operators (still denoted Sg ) acting on functions defined on the boundary OQt o-nd meeting properties (77^- 
|yff[ ), with the norms \\ ■ ||ffQ(asiT)- 

Now, following 13 El], we describe the iteration scheme for problem (Tfi ]) -([76 |) . We choose 

Uo = 0, </?o = 

and assume that {Uk, fk) are already given for fc = 0, . . . , n. Moreover, let (Uk, ^k) vanish in the past, i.e., 
they satisfy ([76]). We define 

where the differences 5Un and Sipn solve the linear problem 

K{U- + U„+i,2. + -^n+mWu = fn in r!r, 

®Ui/2('^^"' ^<^») - 9n on d^T, (80) 

{5Un,5^n) = fori<0. 
Here 

5U^:^5Ur.-——^ rai((/ + C/° + C/„+i/2) (81) 

is the "good unknown" (cf. ([27]) ). 

the operators Lg and are defined in ((29]) . ([30]) . and {Un+i/2,Vn+i/2) is a smooth modified state such 
that {U°- + C/,i+i/2, (yf" + <)5,i+i/2) satisfies constraints (f24]) - (|26]) (^'„, ^'„+i/2, and (55'„ are associated to 
'/'m </'n+i/2i and Sipn like 'J is associated to ^p). The right-hand sides /„ and gn are defined through the 
accumulated errors at the step n. 

The errors of the iteration scheme are defined from the following chains of decompositions: 

£([/„+!, *„+i)-/:(i7„,4'„) 

= L'([/- + C/„+i/2, + *„+i/2)(<5(7„, ,5v|/„) + e'„ + e;^ + e^^' 

= U(C/" + C/„+i/2, ^-"^ + ■^n+l/2)5Un + + + C^^' + i?„+l/2<5^'„ 



and 



=K+i/2isun,spn)+K^+e:+e'::, 
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where Sg^ are smoothing operators enjoying the properties of Proposition l4.1i with the sequence defined 

by 



and we use the notation 



1 



n+l/2 



> 1, 9n = ^/e^+^, 



9l($«+«'„+l/2) 

The errors and are the usual quadratic errors of Newton's method, and e", e" and e"', e"' are the first 
and the second substitution errors respectively. 
Let 



'1 _1_ 



+1/2 



then the accumulated errors at the step n>l are 



En — ^ e-k, En = e-k, 



fe=0 



fc=0 



(82) 



(83) 



with Eq '= and •= 0- The right-hand sides fn and Qn are recursively computed from the equations 

71 n 

; fk + Se„ En = Se„ f , 5] ^fc + Se„ En = 0, (84) 



E. 

fc=0 



fc=0 



where /o :— Segf"' and (70 0. Since Sg^^ — s- / as — > cx), one can show that we formally obtain the 
solution to problem ([7il) - ([75)) from £([/Ar, ^I'jy) ~^ and B{Un\xi=o, ^n) 0, provided that (cat, cat) 0. 

Remark 4.1 In general, the realization of the Nash-Moser procedure for problem (I7^ - ([7S)) below is much 
simpler as in [5T] for current-vortex sheets. As in 5 and unlike [H], we work in usual Sobolev spaces H'' 
(in [5T] one works in the anisotropic weighted Sobolev spaces H^). More precisely, in |5J the exponentially 
weighted Sobolev spaces :— e^^H^ were used, but for H^, Sobolev's embeddings, Moser-type inequalities, 
etc. are internally the same as for the usual Sobolev spaces . Therefore, in some places below our 
calculations are almost the same as in [5]. However, for convenience of the reader we prefer to present all 
the calculations (at least, in brief). Moreover, since, unlike 0, we do not assume that our initial data are 
close to a constant solution and in our tame estimate (|43|) we lose, as [21] , "one derivative from the front" , 
somewhere we have to modify arguments of [S] . 

Below we closely follow the plan of [5] and [2l]- Let us first formulate an inductive hypothesis. As in [21] 
and unlike [5j, we do not require more regularity for 5tpk in our inductive hypothesis. 

Inductive hypothesis. Given a small number (5 > 0, the integer a :— m + 1, and an integer a, our 
inductive hypothesis reads: 

a) Vfc = 0,...,n-1, Vse[3,a]nN, 

¥Uk\\HHnr) + ¥vk\\HHan^) < S9r''-'Ak, 
h) V/e = 0,...,n-l, Vs G [3, a - 2] n N, 



{Hn-l) { 



c) Vfc = 0,...,n-1, Vse[4,a]nN, 

\\B{Uk\.,=^.Vk)\\H'iOnr)<K"'~\ 
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where = 9k+i — Ok- Note that the sequence (A„) is decreasing and tends to zero, and 

Recah that (Uk,ipk) for k = 0, . . . ,n are also assumed to satisfy (|76p . Running a few steps forward, we 
observe that we wih need to use inequalities ((7T|) and ([73]) with m — a — A. That is, we now choose 
a — m + i. Our goal is to prove that (iJ„_i) implies (Hn) for a suitable choice of parameters 9o > I and 
(5 > 0, and for a sufficiently short time T > 0. After that we shall prove (Hq). From now on we assume that 
(Hn-i) holds. As in [5], we have the following consequences of 

Lemma 4.3 If 9o is big enough, then for every k — 0, . . . ,n and for every integer s € [3, a] we have 



\\Uk\\HHnT) + yk\\HHonT)<K' ''>^, a^s, (85) 

\\Uk\\H''{nT) + WfkWH'-idnT) < '^fog^'fc, (86) 

ll(/ - Sg,)Uk\\HHn^) + 11(1 - Sg,)ipk\\Hs(^onr) < CSei-". (87) 
For every fc = 0, . . . , n and for every integer s £ [3, ci + 4] we have 

WSe^UkUnnr) + WSe.^kWH^anr) < CSef'>\ a^s, (88) 

\\Se,Uk\\HHn^) + \\Sg,^k\\H^ian^) < C6\og9k. (89) 



Estimates (|87l)~(189l) follow from and Proposition O Moreover, ([STD and ^ hold actually 

for every integer s > 3 but below we will need them only for s € [3, a] and s € [3, a + 4] respectively. 

Estimate of the quadratic errors. The quadratic errors 

e'fe = CiUk+i,^k+i) - CiUk^'^k) - C'iUk,'^k){6Uk,S-i'k), 
e'fe = {BiUk+i,iPk+i) - B{Uk,iPk) - B'{Uk,^kK6Uk,6ipk)) Ui=o 

can be rewritten as 

e'fe= / {l-T)h"(U'' + Uk + TSUk,^'' + ^k + r5^k){{SUk,S<fk)ASUk,S<fk))dT, (90) 

e'k ^^M"{iSUkU,=o,Sipk)ASUkU,=o,Sipk)) (91) 
by using the second derivatives of the operators L and B: 

V'iU, mU', *'), iU", *")) := :^L'(C/„ ^,){U', ^%=o (L'(C/, $)(C/", *") := -^L(f/„ ^f,)), 

de de 

where Ue — U + eU", We = U\xi=o + sW" , ipe — (p + et/?", and and ^P" are associated to tp' and ip" 
respectively like ^ is associated to <p. We easily compute the explicit form of B", that do not depend on the 
state ([/, 0): 

To estimate the quadratic errors by utilizing representations ([50]) and (PT|) we need estimates for L" and 
B". They can easily be obtained from the explicit forms of L" and B" by applying the Moser-type and 
embedding inequalities. Omitting detailed calculations, we get the following result. 



i"{{W',p'liW",p"))^ 
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Proposition 4.2 Let T > and s eN, with s > 3. Assume that {U,ip) G i/'+i (f^r) x H^'+^idnT) and 

\\U\\H3((nT) + WfllH^idnr) ^ ^■ 

Then there exists a positive constant Kq, that does not depend on s and T , and there exists a constant 
C{Ko) > such that, if K < Kq and {U',(p'), {U",(p") G H^'+^inr) x H^'+^idflT), then 

\\h^uM{U\^')AU'\^n)\\HHnr)<CiKo){((iUj)))^^,m^ 

where (((C/,^))), := ||(7|U'(o,) + Mn'ion^)- If{W',^'), {W",^") G H^idrir) x H^+^dnr), then 

||B''((W-^<^0, < C(i^o){||W^1U^nao.)ll^''IU3(ao.) 

+ \\W\\H3(dnT)\\'P"\\H'+HdnT) + \\W"\\H^(^dnT}\\'P'\\H^dnT) + \\W"\\H^(^dnT)y\\H=+HdnT) 
+ \\W'\\HHdnr)\\W"\\HHan^) + \\W'\\HHdn^)\\W"\\H^^9n^j}. 

Without loss of generality we assume that the constant Kq = 2C*, where C* is the constant from (|7T|) . 
By using ((90)) . ((9T|) . and Proposition [42l we obtain the following result. 

Lemma 4.4 Let a > 4. There exist S > sufficiently small, and 9q > 1 sufficiently large, such that for all 
k = 0, . . .n — 1, and for all integer s G [3, 5 — 1], we have the estimates 

\\ei\\H^^nr^<C5'e^,^^'^-'A,, (92) 
\\e',\\Hsian^^<C5'9t'^'^-'A,, (93) 

where Li(s) = max{(s + 1 — 0;)+ + 4 — 2a, s + 2 — 2a}. 

Proof. In view of (|7ip (recall that m = a — 4), (_ff„_i), and ([85]) . we estimate the "coefficient" of L" in ((90| 
as follows: 

sup ((([/" + iJfc + rSUk, + </?fc + ^'^9'fe))>3 < a + + Sel-^'Ak < C, + C5 < 2C, 

T-e[o,i] 

for S sufficiently small. Therefore, we may apply Proposition [4?2l 

||4|k=(o,) < C[6'9t-''^AI{C. + m,,^,)))^^, + (((<5C/,,<5^,.)L+i) + S'ei+'-'^Al) 
for s G [3, 5 — 1]. If s + 1 ^ a, it follows from that 

(here we have used the inequality 9k Ak < 1/2). If s + 1 = a and a > 4, 

ll4lk=(o,) < C<5^A2 {(a + (51og(?fc + (50-iAfc)^?4-2" + 01-} < C<5^a2(?i-" < CS^^''^''-'^-' A^. 
Analogously, by using ([91]l . Proposition [4?2l and the trace theorem, we get (|93|) . □ 
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Estimate of the first substitution errors. The first substitution errors can be rewritten as follows: 

1 

L" [U'' + Se, Uk + t{I - Se,)Uk, + Se, 'fk (94) 
+t(/ - 5eJ*fc) {{6Uk,6^k), ((/ ~ SgjUk, (/ - Se,)-^k))dT, 
el = {B'{Uk,^k)iSUk,Sipk) - B'{Se,Uk,Se,^k){SUk,Sipk))U^Q 

(95) 

= M"[{6Uk\xi=o,5ipk), {{Uk - S0i^Uk)\xi=o, Vk - Se^ipk))- 

Lemma 4.5 Let a > 4. There exist 5 > sufficiently small, and Oq > 1 sufficiently large, such that for all 
= 0, . . . n — 1, and for all integer s g [6, 5 — 2], one has 

< C5^6l^^^^-'^k. (96) 
\\el\\H^(an^)<C5^e^^^'-^^-'/\k. (97) 

where L2{s) = max{(s + 1 — 0;)+ + 6 — 2a, s + 5 — 2oi\. 

Proof. It follows from jTlI), (-ffn-i), jHT]), and ([88]) that 

sup ((({/" + ^e.L/fe + r(/ - Se,)Uk,v" + Se.^k + t{I - Sejipk)))^ < 2C, 
Te[o,i] 



for 6 sufficiently small, i.e., we may apply Proposition 14.21 for estimating L" in ((94)) . Using again (|7T|). 
(i?n-i), dSZl), and (IMI), for s + 1 7^ a and s + 1 < a we get 

IKllH^in^) < c{s%-''^Ak{c. + de^^+'-''^^ +dei+'-'^)+s%+'-''^Ak} < csHl^^^'^-^Ak. 

Similarly, but exploiting instead of , for the case s + 1 = a we obtain 

l|e'fc'||ff=(n,) < c[5^el-^''Ak{c^+6\ogek + 6) + 5^elr''Ak] 

for a > 4. 

By virtue of ([95]) . the trace theorem, and Proposition [42l we have 

llefellif.(Or) < c[[SUk]s+i,*M\{^ - SeJvkWmidnT) + \\SUk\\m{nT)\\i'^ - 'S'ej¥'fc|U=+i(anT) 
+ ll(^ - Sek)Uk]H-+^nT)\\^'fk\\H^dnT) + ll(^ ^ Se^)Uk\\H^nT)\\^'fk\\H-+^dnT) 
+ \\SUk\\H^+HnT)\\{I - SeJUk\\H^(nT) + ll^t^fcll//3(n^)||(/ - 5ejC/fc||^e+i(f^^)|. 
Then, {Hn-i) and ^ imply □ 



Construction and estimate of the modified state. Since the approximate solution satisfies the strict 
inequalities ([6]) (for all x G fix) and (|2T|) (see Lemma IT2|) and since we shall require that the smooth modified 
state vanishes in the past, the state {U"" + C/„+i/2, + 'Pn+1/2) will satisfy © and (HH) for a sufficiently 
short time T > 0. Therefore, while constructing the modified state we may focus only on constraint (|25p . 
i.e., the first boundary condition in (|17[) . 
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Proposition 4.3 Let a > A. The exist some functions Un+1/2 o^n-d that vanish in the past, and such 

that {W^ + ?7„_|_i/2, + satisfies 125\) . and inequalities 0) and \21\) for a sufficiently short time T . 

Moreover, these functions satisfy 

Vn+1/2 ^ Se^ipn, Pn+1/2 ^ Se^pn, Vj^n+1/2 ^ Se^Vj^n (j = 2,3), Sn+1/2 = Se^S^ (98) 

and 

\\Un+i/2 - Sg,^Un\\H^inr) < 050^+^-" for s e [3, a + 3]. (99) 
for sufficiently small 5 > Q and T > 0, and a sufficiently large Oq > 1. 

Proof. Actually, estimate ([M]) which we are going to prove hold for every s > 3 but below we will need it 
only for s € [3, a + 3]. Let (Pn+i/2j the pressure p„+i/2, the entropy Sn+i/2^ a-nd the tangential components 
of the velocity Vn+1/2 are defined by (fM]) . We define wi^„+i/2 as in |21j : 

t'l,ri+l/2 := Se^Vi^n +'R-tG, 

where 

3 

g = dt(fn+l/2 - {Ser,Vl,n)\xi=0 + ^{{Vj + ^^j,«+l/2)l9j(/3„+l/2 + Vj,n+l/2dj(f'')\^^^Q, 

and TZt ■ H^{d^T) — * H'''^^{VIt) is the lifting operator from the boundary to the interior. To get the 
estimate of wi,„+i/2 — Se^vi^n we use the following decompositions: 

3 

g = Sg„Bi{Un\xi=0,'Pn) - dt{l - Sejfn + (1 ^ <S'e„ )9t (^„+ ^ ((w° + Se^Vj^n)djSg^ipn 



la;i=0 



and 



Bl{Un\xi=0, Vn) = Bu(C^n-lUi=0, V'k-i) + dt{Sipn~-l) 



+ ^((Wj + Vj.7i-l)dj{Sipn-l) + 6Vj^n-ldj{(f'' +(pn) - Svi^n-l) 
3=2 

where Bi denotes the first row of the boundary operator B in (I75p . 
Exploiting point c) of (_ff„_i), one has 

\\TiT{S0^Bi{Un~l\xi=O, Vn-l))\\H=(nT) < C\\Se„Bi{Un-l\xi=0, ^n-l)\\H=(dnT) 

J C6'f-"||6i(f7„_iUi=o,'^«-i)|!H=(anT) for s e [a,a + 3], , oi 

<< < C(56'^ " for s e [3, a + 3]. 

y C\\Bi{Un-i\xi=o,Vn^i)\\H-+\a^iT) for s G [3, a - 1] 

Using (|77| and point a) of (_ff„_i), we get 

pi-T{Se,A{^V>n-i))\\H- {Ut) < C\\Se^dt{5<~Pn-i)\\H-(dnT) 

for s S [3, a + 3]. We also obtain 

W^riSeAiv^ + «j,n-i)\x,=odjiSif,,^i)))\\H^^nT) < C^r'IK^" + ^'j,«-i)Ui=o a,((5(^„_i)|U3(ao^) 
< Cdfr''{\\SiPn-i\\mionr)\\U'' + C/„-i||h3(o,) 
+ \\S^n-i\\HHanr)\\U'' + U^-i\\H^(nr)} < COf^-'Sel-^'C, < C5e--<-^ 
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for j = 2, 3 and s G [3, a+3]. Estimating similarly the remaining terms containing in TZTiSo^BilUnlxi^o, 
we finally obtain 

\\nT{Se„Bi{UnU,=o,VnmHHnT) < CSe:-", se [3,5 + 3]. 
We now need to derive estimates for the remaining terms containing in TZtG- For s G [a, a + 3] one has 

|j7^T(-l9t(l - SeJ(pn + (1 - S0jdt(pn)\\H'(nT) < C{\\dt{Se^(pn)\\H-(anT) + WSe^idt^Pn^H-ianT)} 
while for s G [3, a — 1] we obtain (recall that a — a + 3) 

||7^T((1 - SeJdtVn)\\H^in^) < C0r"||<^„|U.+i(aa,) < CJ^'""- 



Here we have, in particular, used Lemma l4.3l We do not get estimates for all the remaining terms containing 
in TZtG and leave corresponding calculations to the reader. Collecting these estimates and the estimates 
above, we finally have 

l|wi,n+l/2 - <5'0„t'l,n||ff=(J2T) < C(56',^+^"" , SG [3, Of + 3], 

that is equivalent to (|99p . □ 



Estimate of the second substitution errors. The second substitution errors 
and 

eJc" = {B'{Se^Uk,Se^^ipk){SUk,Sifk) - l3\Uk+i/2,(Pk+i/2){^Uk,Sipk))\xi=o 

can be written as 

e'fc"= / V'iU'' + Uk+i/2 + r{Se,Uk-Uk+i/2),^'' + Se,'i'k){{SUk,S-^k),{Se,Uk-Uk+i/2,0))dT, (100) 
Jo 

4' =M"{{SUkU=o,S^k)AiSe,Uk-Uk+i/2)\x^=Q,0))- (101) 
Employing (|100p and (jlOip . we get the following result. 

Lemma 4.6 Let a > A. There exist 6 > 0, T > sufficiently small, and 9q > 1 sufficiently large, such that 
for all k — 0, . . .n ^ I, and for all integer s G [3, 3 — 1], one has 

\K\\H^(^n.^<CS'e;:^^^^-'Ak (102) 

and e'lf ~ 0, where Lz[s) — max{(s + 1 — «)+ + 8 — 2a, s + 5 — 2a}. 

Proof. Using Lemma 14.31 and Proposition 14. 3|, we obtain the estimate 

sup i(U<' + + T(Sg,Uk - Uk+1/2), + Se.^k)))^ < 2C, 

re[o,i] 
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for 6 sufficiently small, i.e., we may apply Proposition [4?2l Similarly, one gets 
((([/" + C/fe+i/2 + T{Se,Uk - Uk+i/2),V^ + Se^^k)))^^, 

Applying Proposition [121 we obtain (|102p : 

Using the explicit form of B", we easily get e'^" = 0. □ 
Estimate of the last error term. We now estimate the last error term 

+ ^'„+i/2) 

where Rk di {L([/" + C/fe+i/2,^'" + *fe+i/2)}- Note that 

+^r.+ l/2)\ = |l+ai(*'^ + *„+i/2)| > 1/2, 

provided that T and S are small enough. 

Lemma 4.7 Let a > 5. There exist 6 > 0, T > sufficiently small, and 9q > I sufficiently large, such that 
for all k — 0, . . .n — 1, and for all integer s G [3, 3 — 2], one has 

\\Dk+i/2S^k\\H^in^) < CS'et^'^^'Ak, (103) 
where L{s) = max{(s + 2 - a)+ + 8 - 2a, (s + 1 - a)+ + 9 - 2a, s + 6 - 2a}. 

Proof. The proof follows from the arguments as in [HIS] (see also [3T])- Using the Moser-type and embedding 
inequalities, we obtain 

\\[Dk+i/2S^k\\H-{nT) ^ C!\ \\6ipk\\H-{dnT)\\Rk\\H3{nT) 

^ (104) 

+ \\S'Pk\\HHanT){\\^k\\H^{nT) + \\Rk\\HHnT)\\f'' + fk+i/2\\HHdnT)) J 

(note that ||9i(^'"+^'„+i/2)||h=(Ot) — ^\\v^ + 'fk+i/2\\H={dnT))- To estimate i?^ we utihze the decomposition 

L(C/- + Uk+i/2, *° + ^k+1/2) = CiUk, ^k) ~r+ L(C/° + C/fe+1/2, + *fe+i/2) 

-L(C/" + Uk, *° + *fc) = C{Uk, *fc) - r+ / L'(C/'' + Uk+ r{Uk+i/2 - Uk), 

Jo 

+ ^-fc + T{^k+l/2 - */c)) iUk+1/2 - «'fc+l/2 - '^k)dT. 

Clearly, 

\\R\\Hs(nr)<\mUk,^k)~n\H^inr)+ sup ||L'(...)(...)||ff.+i(n,) (105) 

Te[o,i] 

(for short we drop the arguments of L'). It follows from point b) of (iJ„_i) that 

\\£iUk, ^k) - r ||ff=+i(n,) < 2<50r" (106) 
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for s E [3, a — 3]. We estimate L' similarly to L" (see Proposition 14. 2[) . One has 

sup ((([/" + Uk+ T{Uk+i/2 - Uk), + + T{^k+i/2 - Vk)))). < 2C, 
■re[o,i] 

for S small enough. Then, omitting detailed calculations, we get the estimate 

||L'(. ..)(.. .)\\Hs+Hn^) < CS{ei+'-" + 

for s e [3, a - 3]. This estimate, (ITU5)) . and (fTITO)) imply 

\\R\\H^in.) < CS{el+'-'^ + (107) 

for s E [3, a — 3]. For s = a — 2 we estimate as follows: 

\\R\\HHnT) < WHU"" + Uk+1/2, + *fc+i/2)||H=+i(n^) 

< C(((C/« + {Uk+i/2 - Se^Uk) + Sg^Uk, f'^ + Se^Vk)))^^^ < 
That is, we get estimate (|107p for s E [3, a — 2]. Using then (|104p . we obtain (|103p . provided that a > 5. □ 

Convergence of the iteration scheme. Lemmas 14 .4114 . 71 yield the estimate of e„ and e„ defined in ([5^ 
as the sum of all the errors of the fcth step. 

Lemma 4.8 Let a > 5. There exist 6 > 0, T > sufficiently small, and 9o > 1 sufficiently large, such that 
for all k — 0^ . . .n ~ 1, and for all integer s E [3, 3 — 2], one has 

WekWHHnT) + WekWHHdnT) < CS^Oj:^'^'' Ak, (108) 
where L(s) is defined in Lemma \4- ?[ 

Remark 4.2 In principle, we could try to use the advantage of the fact that in the tame estimate P5|) we 
do not lose derivatives from the source term / to the solution. To this end, in Lemma 14.81 we could estimate 
errors e„ and e„ separately. However, this does not reduce the number of derivatives lost from the initial 
data to the solution in the existence Theorem 1 2. II In fact, we can even use a roughened version of estimate 
(|43p in which we lose one derivative from / to the solution. 



Lemma [4.81 gives the estimate of the accumulated errors En and En. 
Lemma 4.9 Let a > 7. There exist S > 0, T > sufficiently small, and 9o > 1 sufficiently large, such that 

ll^'„||//a+2(o^) + \\En\\H''+-^{dnT) ^ CS'^On, (109) 
where L{s) is defined in Lemma \4- ?[ 

Proof. One can check that L{a + 2) < 1 if a > 7. It follows from (fTU5|) that 

n—1 n—1 

{{{En,En)))^+, < (((efc,efe)»„+2 < E ^ 

A;=0 fe=0 

for a > 7 and a + 2 E [3, a — 2], i.e., a > a + 4. The minimal possible a is a + 4, i.e., our choice a = a + A 
is suitable. □ 

We now derive the estimates of the source terms /„ and (?„ defined in ([84]) . 



27 



Lemma 4.10 Let a > 7. There exist 6 > 0, T > sufficiently small, and Oq > I sufficiently large, such 
that for all integer s G [3,3 + 1], one has 

WfnWHHnr) < CA„{0r"-2 (||r||«„+i(^,,) + 6') + (110) 

hnWH^ian^) < CS'A^e^^^^-' + e^-^-'). (Ill) 



Proof. It follows from ^ that 

fn — (Se,, — 'S'e„_i)/" — {Se^ — Se^_i)En-i — S'e„e„_i. 
Using ([77]), dlSl), (UnS), and pUg]) . we obtain the estimates 

||(5'e„ - S0^_^)En-l\\\H= (ilT) < C'^T^iIl~^l|-E'n-l||//°+2(i2T)^n-l - C''^^ ^n- ^n- 1 : 

Using the inequalities 6'„_i < On < \/2^?i-i, ^^n-i < 30„, and A„_i < 3A„, from the above estimates we 
deduce pTU)) . Similarly, we get pTT]) . □ 

We are now in a position to obtain the estimate of the solution to problem (|80p by employing the tame 
estimate P5)) . Then the estimate of {SUn,Sipn) follows from formula (|ST|) . 

Lemma 4.11 Let a > 7. There exist 6 > 0, T > sufficiently small, and Oq > 1 sufficiently large, such 
that for all integer s € [3, a], one has 

\\SUjHs(n^) + \\Sipn\\Hs(anr) < <5^?r"-iA„. (112) 



Proof. Without loss of generality we can take the constant Kq appearing in estimate (|43p that Ko = 2C* 
where C* is the constant from (|7T|) . In order to apply Theorem 13. li by using (|88p and ([99)) . we check that 

IIC^" + C^ri+i/2llH6(j2T) + ll'y'" + '5'e„'Pri||fl-6(af2T) ^ 2C, 



for a > 7 and (5 small enough. That is, assumption (|42)) is satisfied for the coefficients of problem ((80)) . By 
applying the tame estimate (|43ll . for T small enough one has 



\\SUn\\H={nT) + W^VnllH^idnT) ^ C'|ll/n||-ff»(J2T) + 1 1 3" 1 1 -ff + i (90t) 

+ (ll/ri||H3(f2T) + \\9n\\HHdnT)){\\U'' + 1 1 //=+3 (J2t) + 11'/'° + Se^^PnW H-+^dnT))} ■ 

Using Moser-type inequalities, from formula (|8ip we obtain 

\\SUn\\H-{nT) ^ ll^f^n||H=(f2T) + C"! 1 1 '^'/'n 1 1 » (SOt ) + 1 1 ^V'" II //^ (SOr) II + 'S'e„ V'n llff= (a^T) } • 

Then pl^ yields 

||'5?7,i||//=^(f2r) + W^^nWH-idnr) < II /n || H= (Or) + \\9n\\ H= + ^ {dnr) 

+ (ll/n||H3(Or) + ll.9n||ff''(ai2T)) (l|t^°' + 1 1 H=+3 (Ht) + II'/'" + '5'e„ Vn 1 1 + 3 (ai2T) ) } 



(113) 



(114) 
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for all integer s e [6, a]. Below we can actually use a roughened version of (|114[) (see Remark 14.21) . Applying 
Lemma [4. Ill (|88|) . and Proposition [43l from (|114p we derive the estimate 

(115) 

Exactly as in [5 , , we can check that the inequalities 

L{s + I) < s — a, (s + 3 — a)+ + 2 — a<s — a — 1, 

(s + 3 - a)+ +9 - 2q! < s - a- 1, (116) 
s + 6 — 2a < s ~ a ~ 1, s + 13 — 3a<s — a — 1 

hold for a > 7 and s e [3, a]. Thus, PT5)) and ^ yield 

P!7„||Hno.) + ll'5^n||H-^(ao,) < C {So{T) + S') ^r^-^A,, < SO^-^-'A^ 
for (5 and T small enough. □ 

Remark 4.3 As we can see, Lemma [4. Ill with a = a + 4 is absolutely analogous to Lemma 16 in [5]. In this 
sense, the "gain of one derivative for the front" in the tame estimate gives no advantage in the realization of 
the Nash-Moser method. This is caused by the fact that even if in point a) of {Hn-i) we had the iJ'^+^-norm 
of Sipk we could never use this advantage before the proof of Lemma 14.111 

Inequality (|112p is point a) of (i/„). It remains to prove points h) and c) of (iJ„). 

Lemma 4.12 Let a > 7. There exist 6 > 0, T > sufficiently small, and Oq > I sufficiently large, such 
that for all integer s G [3,5—2] 

\\C{Un,^n)~flHHnT)<2S0:,~"~'- (117) 
Moreover, for all integer s G [4, a] one has 

\\B{UnU,^o, ^n)\\HHan^) < Set'^'K (118) 



Proof. One can show that 

C{Ur,, *„) - r - {Se,^_, - I)r + (/ - 5e„_ + e„-i. (119) 

For s G [a + 1, a — 2], by using ([77]) . we obtain 

W~Se^_,)f^\\H^(n^) < ^r"-i(C||r|U.+i(a,) + ||rik.(a.)) < C5o(r)^?r"-\ 
while for s G [3, a + 1], applying ([78|) . we get 

\\{I-Se„^,)f^HHnr) < C0:-r'\\n\H^^Hnr) < C5o{T)e:-"-\ 
Lemma SSI and ([75)1 imply 

Wil - Se„_jE„-l\\H^{nT) < C'^n-l ^l|-E'ri-l||_H-" + 2(Or) < CS'^6^^°'~^ 
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for3<s<a + 2 = a- 2It follows from (fT08)) that 

From the above estimates and decomposition (|119p . by choosing T > and 5 > sufficiently small, we 
obtain pi7p . Similarly, by using the decomposition 

B{Un\xi=0, fn) = U - Se^_i)En-l + e„_i, 

we can prove estimate (jllSp . □ 

As follows from Lemmas 14.111 and 14 . 1 2( we have proved that implies (Hn), provided that a>7, 

a = a + A, the constant 0o > 1 is large enough, and T > 0, S > are small enough. Fixing now the constants 
a, 5, and Oq, we prove (Hq). 

Lemma 4.13 If the time T > is sufficiently small, then (Hq) is true. 

Proof. We recall that {Uojfo) = 0. Then, by the definition of the approximate solution in Lemma [4.21 the 
state ([/" + Uq, (f°- + (fa) — satisfies already (fT7| . and ([2T|l . That is, it follows from the construction of 
Proposition 14. 31 that {U1/2, — 0. Consequently, ((5{7o, Sipo) solves the linear problem ([3T |) - ([33|) with the 
coefficients {U,(p) — {U'^,(p°') and the source terms / — Se^f"" and g — 0. Thanks to (|7ip the assumption 
(|42| is satisfied (recall that Kq — 2C*). Applying (|43p . we get the estimate 

ll'5C/o||_fi-=(OT) + \\^Vo\\H=(dnT) ^ CWSegf^Wn'+^QT)- 
Together with ([74|) and formula (l8T|) this estimate yields 

\\SU,\\H^,(nr) + WS^H^idnr) < C||5eor||^f.+i(o,) < C0^^""^+(5o(T) < 69',-"-' Ao 

for all integer s G [3,(5], provided that T is sufficiently small. Likewise, points b) and c) of (Hq) can be 
shown to be satisfied for a sufficiently short time T > 0. □ □ 

The proof of Theorem HHl We consider initial data {Uo,(po) G x H™+'^(R2) satisfying 

all the assumptions of Theorem 12.11 In particular, they satisfy the compatibility conditions up to order 
~ m + 7 (see Deffiiition 14. ip . Then, thanks to Lemmas 14.11 and 14.21 we can construct an approximate 
solution iU'',ip'') e H"^+^{flT) X H"^+^{dnT) that satisfies As follows from Lemmas IHlHlTl (7J„) 

holds for all integer n > 0, provided that a>7,a — a + A, the constant 6*0 > 1 is large enough, and the 
time T > and the constant 6 > are small enough. In particular, (i/„) implies 

00 

{WSUnWH^iflT) + \\S^n\\H"^idnT)} < 

Hence, the sequence ({/„, (fin) converges in H"^{VLt) x H^^{dilT) to some limit ([/, (p). Recall that m = a — 1 > 
6. Passing to the limit in (fTTT]) and (fTTS]) with s = m, we obtain (I74p - (l76p . Consequently, U := U + U"', 
(f :— ip + (f^ is a. solution of problem (^0]) . pT]) . This completes the proof of Theorem 12. II 
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5 Free boundary problem in relativistic gas dynamics: special and 
general relativity 

Let us first write down a suitable symmetric form of the relativistic Euler equations. First of all, we note 
that for the set of covariant laws ([5]) we have the supplementary covariant law 



V„(pS'u") = 



(120) 



that arises as a consequence of ([8]) and the first principle of thermodynamics. In the setting of special 
relativity l|120p becomes the entropy conservation law 



9t(prS') +div {pSu) = 0. 



(121) 



In principle, taking into account p2ip and using Godunov's symmetrization method, we can rewrite system 
(|10p - p2p for the unknown U — (p, u, S) as a symmetric system for a new (canonical) unknown Q and then 
return to the original unknown U keeping the symmetry property: 



A"{U)dtU + A^iU)djU + Q{U) = 0, 



(122) 



where A" = (A")^, dj — d/dx^, and Q{U) — — (0, — pC/, 0). This procedure is described in [3] where the 
symmetric matrices A" were written for the special case = u"^ = 0. Such a procedure is absolutely 
algorithmic and always works, but it is however connected with very long calculations. Therefore, here we 
prefer to symmetrize the conservation laws pII)) - ([T^ by rewriting them in a suitable nonconservative form. 
Equations ^ and p2T|) imply 

f.O, ,123) 

where d/dt — dt + {v, V) is the material derivative as for the non-rclativistic case (|4]). Using (|123p . we first 
rewrite pU]) in a nonconservative form. Combining then pT|) and ([T^) and employing again (|123p . we finally 
get the relativistic counterpart of system (|3|): 



% + {v, dtu) + div u = 0, 
p& at 



{phT) 

dS_ 
dt 



du 
dt 



dt 



{dtp)v + Vp = pQ, 



(124) 



0, 



where c — {pp{p, S)) 



1/2 



System (jl24p being written in the quasilinear form (|122p is already symmetric with 



pc- 

V 



phTL 




A^ = 



( 









\ 




pc^ 










phu^Sg 







\ 










) 



(125) 



where 3§ = (bij), hi 



{Sij,S2j,S3j), and is the vector-row for a corresponding 



vector-column a (recall also that = Tv^). The matrix Aq > provided that inequalities ^ are satisfied 
together with the relativistic causality condition 



< c: < 1, 



(126) 
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where Cg is the relativistic speed of sound, = c^/h. Of course, (|126p will be an additional restriction on 
the initial data in a counterpart of Theorem 1 2. II 

Now, for system (|122|) . ()125|) in the domain (fT3|) endowed with the boundary conditions ([T4|) we can 
literally repeat arguments of Sections UHll The only important point is that the boundary matrix Ai on the 
boundary xi = for system (|34)) written now for matrices (|125p and V — {pjUmU^ , S) coincides with 
the matrix y^i|a;i=o in (|35p . where 



Un-.^fvn, Vn=i^ -V^d2^ -V^ds^, f = (l + \u\Y/^ , V=u/f, 



(127) 



U — {p, ii, S) is the "good unknown" , U — {p, ii, S) is the basic state, and v = (v^ , w^, v^) is defined from the 
formula 

u = fv + fu{u,v) (128) 

suggested by the relation between the perturbations Su and Sv. 
Indeed, we easily compute: 







a phT^^ 
V f / 



where 



a = (1 - vi^t^, -92$ - W2at$, -93$ - i)39t$), f = - w^92$ - w^Sa* - 9*$, 

and ^ is the matrix ^ calculated for the basic state. Taking into account ([25]), (|127[) . and ((l3 

f|a:i=o = and 

(5i$) (Ii(^, *)f7, f/)U,=o = 2 pU^^o (^ti - u''d2^ - 7i'93$ - (i), u)dt^) 



we have 



a;i=0 



2(pWn)Ui=0 = (^(l)V^Ui=0,"l^Ui=o) 



= 2 pU,=o + (u, «)at^(r2 - 1 - |{ip)) 

(the matrix A^i^ was defined in (|35p ). Then 

(Ii (&,$)[/,[/) U,=o = (Ii(&,$)Ji^, Ji^)Ui=o = (J^Ii(&,$)jy,y)Ui=o = (AV^,F)U,=o, 

where the matrix vAil^^j^o is the same as in Section [5] and the transition matrix J can be easily written 
down. Thus, we obtain the local-in-time existence (and uniqueness) theorem for the relativistic version of 
problem pBjl - p^ (in the framework of special relativity) in the form of Theorem 12.11 Clearly, we should 
also supplement conditions ^ with p26p while writing assumptions on the initial data. It means that the 
initial data should satisfy 

inf {p{po,Sa), pp{po, So), c^(po, •S'o), 1 - c^(po, <S'o)} > 0, 



where 



c^(Po,S'o) 



1 



Kpo, So) = 1 + e{p{po, So), Sq) 



Po 



Pp{po,So)h{po,So)' -"v-f"'""^ - ' "vrv^-ui-uy,-.,/ i pi^p^^s^y 
Let us now briefly discuss the case of general relativity. The metric g appearing in the relativistic Euler 
equations ([9]) should satisfy the Einstein equations Ga/3 = uTap. As in [6], following Rendall [TQ and 
introducing 
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we write the Einstein equations in harmonic coordinates as 

g'^dtgapi - g^^digai3z = 0, (129) 
dtgaf3 — 5q/30 = 0. 

System (|129p written in the compact form 

B°(W)dtW + B^{W)djW + Q{W, U) = (130) 

is symmetric for the vector W whose components are gaf3 and gap-y- Recall that U — {p, u, S). The symmetric 
system (|130p is hyperbolic if 5"° < and (5*^ ) > 0. 

Regarding the relativistic Euler equations ([9]), it is enough to symmetrize them for a fixed constant 
metric g. This was done by Kendall |18j for isentropic fluids. In the general case we can however just repeat 
arguments from [18] by taking into account the entropy law (|12ip which has form (|123p for constant metrics. 
Roughly speaking, the calculations in [T^ are just a "tensor" variant of our simple calculations towards 
obtaining the nonconservative form (|124p . With reference to [18] , we write equations ([9]) for a fixed constant 
metric g in the symmetric form p22p . p25p with 

^ = (^ii): ^ij = gi] + goiV^ + gojv' + goav'v^ = gij + goi— + goj— + goa-r~oT2- 

For a non-fixed metric g the balance laws ([9]) are written as the symmetric system 

A^{U)dtU + A^{U)djU + B{U, W) = 0. (131) 

It is worth noting that for system (|13ip for any fixed (and not necessarily constant) metric we can prove a 
counterpart of Theorem 12. II under suitable assumptions on W . 

Now we consider the free boundary problem for the symmetric hyperbolic system (|131l) . (jl30p with the 
boundary conditions ()14p . However, in the setting of general relativity it is actually an interface problem 
because we should consider system ()130p for the metric variables not only in the domain il(i) but also in 
the vacuum region M^\51(t) — {x^ < ip{t, a;^, x^)}. As was shown in [8 , the jump conditions on an interface 
I](i) written for the Einstein tensor are satisfied if the metric g is smooth on this interface, i.e, 

[W] = W+ -W- ^0 on E(i). (132) 

In our case W'^ and W" are the metric variables in the fiuid domain n(t) and the vacuum region R'^\f2(t) 
respectively. Constraints on the initial data under which condition (|132p is not only sufficient but also 
necessary for the fulfillment of the jump conditions for the Einstein tensor are discussed in [6] and connected 
with the notion of so-called natural coordinates [8]. That is, as for shock waves in general relativity studied 
in [6], we will treat our problem in harmonic natural coordinates. 
Thus, we have the symmetric hyperbolic systems 

A"{U)dtU + A^{U)djU + B{U,W+) = Q mn{t), (133) 
B"{W+)dtW+ + B^W+)djW+ + Q{W+,U) ^0 in f7(t), (134) 
B°{W-)dtW- + B^W+)djW' +Q{W-,0) ^0 mR^\n{t) (135) 

endowed with the boundary conditions (fTH) and (|132p on a time-like hypersurface 'S{t) — {x^ — (p{t, x'^,x^)}. 
Here (|135p is the symmetric form of the vacuum Einstein equations. We reduce problem (|133p - (|135p . (O, 
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(|132p to the fixed domain M.^ by straightening the free surface E: 

U{t,x) U(t,<^>+{t,x),x'), := W^{t,<i>^{t,x),x') 

±a;^ + "ii^it, x) := xi±xi)ip{t, x'), x'^{x'^,x'^) 

(the cut-off function xi^i) was described in the beginning of Section[21). 

Regarding further arguments towards the proof of the local-in-time existence theorem for the reduced 
problem in the domain M5[,, we give here only a rough scheme or even an idea of this proof and postpone 
detailed arguments to a future work. The main idea is the following. The existence of solutions of problem 
(|133p . (HH) reduced to the fixed domain M.^ is proved by Nash-Moser iterations for any fixed metric g. The 
boundary conditions l|132[) are linear and, therefore, we do not need introduce source terms for them in the 
linearized problem. Moreover, for the linearized problem these boundary conditions are dissipative. Though, 
they are not strictly dissipative, but the crucial point is that they are homogeneous. Hence, we can prove the 
existence of solutions to the reduced problem for (|134p . (|135p . (|132p in [0, T] x by the classical fixed-point 
argument for any fixed fluid unknown U . Then, the existence of solutions to the whole problem l|133p - (|135p . 
(HH) , (|132p reduced to the fixed domain M.^ is proved by Nash-Moser iterations for the "fluid" part of the 
problem whereas at each Nash-Moser iteration step the metric g is found as a solution of the problem whose 
linear version has maximally dissipative boundary conditions. More presicely, at each (n + l)th iteration 
step before solving the linear problem for SUn with = W^f we find as a unique solution of the 
corresponding problem for with U ~ Un and ip — (pn taken from the nth iteration step. At last, we 
note that the constraints on the initial data connected with the introduction of natural coordinates are 
not needed to be satisfied at each Nash-Moser iteration step and we may therefore not care about them. 
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